


use map.dta,clear
merge 1:1 center using itdb_prov

spmap years_occupation using itcoord_prov, id(center) fcolor(Blues) clmethod(u) clbreaks(0  2 3 4 5 6 7 8 9) label( xcoord( x_stub) ycoord(y_stub)  label( VARNAME_1 ) size(*0.5 ..) position(0 6) length(26))



